Evidence for long-term Gamma-ray and X-ray variability from 
the unidentified TeV source HESS J0632+057 



V. A. Acciari 1 , E. Aliu 2 , T. Arlen 3 , M. Beilicke 4 , W. Benbow 1 , D. Boltucb 2 , S. M. 
Bradbury 5 , J. H. Buckley 4 , V. Bugaev 4 , K. Byrum 6 , A. Cannon 7 , A. Cesarini 8 , A. 
Cesarini 8 , Y. C. Chow 3 , L. Ciupik 9 , P. Cogan 10 , R. Dickherber 4 , C. Duke 11 , T. Ergin 12 , A. 
Falcone 13 , S. J. Fegan 3 , J. P. Finley 14 , G. Finnegan 15 , P. Fortin 16 , L. Fortson 9 , A. Furniss 17 , 
K.Gibbs 1 , G. H. Gillanders 8 , J. Grube 7 , R. Guenette 10 , G. Gyuk 9 , D. Hanna 10 , J. Holder 2 , 

D. Horan 18 , C. M. Hui 15 , T. B. Humensky 19 , A. Imran 20 , P. Kaaret 21 , N. Karlsson 9 , M. 
Kertzman 22 , D. Kieda 15 , J. Kildea 1 , A. Konopelko 23 , H. Krawczynski 4 , F. Krennrich 20 , M. 
J. Lang 8 , S. LeBohec 15 , G. Maier 10 '*, A. McCann 10 , M. McCutcheon 10 , J. Millis 24 , P. 
Moriarty 25 , R. Mukherjee 16 , R. A. Ong 3 , A. N. Otte 17 , D. Pandel 21 , J. S. Perkins 1 , D. 
Petry 26 , M. Pohl 20 , J. Quinn 7 , K. Ragan 10 , L. C. Reyes 27 , P. T. Reynolds 28 , H. J. Rose 5 , M. 
Schroedter 20 , G. H. Sembroski 14 , A. W. Smith 6 , D. Steele 9 , M. Theiling 1 , J. A. Toner 8 , A. 
Varlotta 14 , V. V. Vassiliev 3 , S. Vincent 15 , R. G. Wagner 6 , S. P. Wakely 19 , J. E. Ward 7 , T. 
C. Weekes 1 , A. Weinstein 3 , T. Weisgarber 19 , D. A. Williams 17 , S. Wissel 19 , M. Wood 3 



ABSTRACT 

HESS J0632+057 is one of only two unidentified very-high-energy gamma-ray sources which 
appear to be point-like within experimental resolution. It is possibly associated with the massive 
Be star MWC 148 and has been suggested to resemble known TeV binary systems like LS I +61 
303 or LS 5039. HESS J0632+057 was observed by VERITAS for 31 hours in 2006, 2008 and 
2009. During these observations, no significant signal in gamma rays with energies above 1 TeV 
was detected from the direction of HESS J0632+057. A flux upper limit corresponding to 1.1% 
of the flux of the Crab Nebula has been derived from the VERITAS data. The non-detection by 
VERITAS excludes with a probability of 99.993% that HESS J0632+057 is a steady gamma-ray 
emitter. Contemporaneous X-ray observations with Swift XRT reveal a factor of 1.8 ± 0.4 higher 
flux in the 1-10 keV range than earlier X-ray observations of HESS J0632+057. The variability in 
the gamma-ray and X-ray fluxes supports interpretation of the object as a gamma-ray emitting 
binary. 

Subject headings: acceleration of particles binaries: general gamma rays: observations stars: individual 
( HESS J0632+057, MWC 148 ) 
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1. Introduction 

HESS J0632+057 is one of about 20 very-high- 
energy (VHE) gamma-ray sources with no known 
counterparts at other wavel engths (for a recent 
review see e.g. IWeeked (120081 )1. Gamma-ray emis- 
sion was discovered by the High Energy Stereo- 
scopic System (H.E.S.S.) during observations of 
the Monoc eros Loop Supernova R emnant in 2004 
and 2005 (jAharonian et all l2007h . It appears to 
be point-like within experimental resolution; the 
limit on the size of the emission region was given 
as 2' (95% confidence level). The reported flux 
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of gamma rays with energies above 1 TeV from 
HESS J0632+057 corresponds to about 3% of the 
flux of the Crab Nebula, with a differential pho- 
ton spectrum consistent with a power-law function 
with index of 2.53 ± 0.26 s t a t =1= 0.20 svs . Possible as- 



sociations considered by lAharonian et al ( 2007 ) 
are the Monoceros Loop Supernova remnant, the 
weak X-ray source 1RXS J063258.3+054857, the 
BOpe-star MWC 148 (HD 259440) and the uniden- 
tified GeV gamma-ra y source 3EG J0634+0521 



(jHartman et al 
servations with 



1999). 



Follow-up X-ray ob- 
XMM-Newton by iHinton et al 
(2009) revealed a variable X-ray source (XMMU 
J063259.3+054801) with a position compatible 
with HESS J0632+057 and MWC 148. It should 
be noted that 3EG J0634+0521 is absent in the 
EGR catalogue of EGRET gamma-ray sources, a 
reanalysis of the EGRET data with new Galactic 
interstellar emission models and interstellar radia- 



tion field data (jCasandiian fc Grenierll2008l ) . 3EG 
J0634+0521 is, as expected from the reported 
flux, not in the Ferm i bright gamma-ray source 
list (|Abdo et alll2009f ). 

Point-like gamma-ray sources stand out among 
the many galactic VHE objects with spatially ex- 
tended gamma-ray emission. The latter are usu- 
ally associated with either pulsar wind nebulae 
or supernova remnants. High-mass X-ray bina- 
ries constitute the only known class of galactic 
objects with variable point-like VHE emission; 
this class currently cont ains three members only : 
PSR B19 59-63/SS 2883 (lAharonian et all l2005bh. 
LS 5 039 (lAharonian et all l2005cl) and LS I +61 
303 (|Albert et all 120061 : lAcciari et all 120081 ) . Ad- 
ditionally, marginal evidence at a level of 3.2cr 
from the bl ack-hole binary Cy g X-l has been 
reported by lAlbert et a3 (|2007fh TeV binaries 
show variable emission of gamma rays, likely con- 
nected to changes in physical parameters associ- 
ated with the orbital movement. VHE gamma- 
ray production in these objects is explained by 
the acceleration of charg ed particles in accretion- 
powered relativistic jets (TTavlor fe Gregory 11984 
Mirabel fc Rodrigue 3 Il994f ) or in shocks cre- 



ated by the collision of the expanding pulsar 
wind with the wind from the stellar companion 
( Maraschi fc Treves] 1981 ). Subsequent inverse- 
Compton scattering on low-energy stellar photons 
produces gamma rays at GeV and TeV energies. 
While there has been no compact companion dis- 
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covered for MWC 148, the point-like nature of the 
VHE emission combined with the variable X-ray 
emission can easily be explained by a production 
scenario similar to those in TeV binaries. A second 
possible scenario is that MWC 148 is a representa- 
tive of a new type of VHE e mitter as proposed by 
Babel & Montmerld (|l997l ) and llbwnsend et all 
(2007h. In their picture strong magnetic fields 
around the massive star lead to magnetically chan- 
neled wind shocks in which second-order Fermi 
acceleration might occur. However it is not clear 
if the circumstellar environment of MWC 148 is 
strongly magnetized, or if this acceleration mech- 
anism is able to produce particles of sufficiently 
high energy to produce a measurable TeV flux. An 
association of HESS J0632+057 with the Mono- 
ceros loop SNR is unlikely given the point-like 
nature of the gamma-ray emission and the non- 
correlation of possible tar get material with the p o- 
sition of the VHE source ( Aharonian et alll2007 ). 



2. Observations 

VERITAS is an array of four imaging atmospheric- 
Cherenkov telescopes located at the Fred Lawrence 
Whipple Observatory in southern Arizona. It 
combines a large effective area (up to 10 5 m 2 ) 
over a wide energy range (100 GeV to 30 TeV) 
with good energy (15-20%) and angular (« 0.1°) 
resolution. The field of view of the VERITAS 
telescopes is 3.5°. The high sensitivity of VERI- 
TAS enables the detection of sources with a flux 
of 1% of the Crab Nebula in less than 50 hours of 
observations. For more details on the VERITAS 
instrument, see e.g. lAcciari et a 3 (|2008h . 

VERITAS observed the sky around HESS 
J0632+057 during three periods in December 
2006, December 2008 and January 2009; see Ta- 
ble [1] for details. For each period, data equivalent 
to 10 hours of observations passed quality selec- 
tion criteria, which remove data taken during bad 
weather or with hardware-related problems. Data 
were taken on moonless nights in wobble mode, 
wherein the source was positioned 0.5 degrees from 
the camera center with the offsets in different po- 
sitions for different runs. The first data set (Set 
I) consists of observations taken during the con- 
struction phase of VERITAS with only 3 tele- 
scopes. These observations were pointed towards 
the centre of the Monoceros region (at an an- 



gular distance of ~ 0.5° from HESS J0632+057), 
while observations in the second and third set were 
targeted around the reported position of HESS 
J0632+057. 

The data analysis steps consist of image cali- 
bration and cleaning, s econd-mome nt parameter- 
ization of these images ( Hillasl 1985). stereoscopic 
reconstruction of the event impact position and 
direction, gamma-hadron sepa ration, spectral en- 



ergy reconstruction (see e.g. iKrawczvnski et al 
( 20061 )) as well as the generation of photon sky 
maps. The majority of the far more numerous 
background events are rejected by comparing the 
shape of the event images in each telescope with 
the expected shapes of gamma-ray showers mod- 
eled by Monte Carlo simulations. These mean- 
reduced- scaled width and mean- reduced- scaled 



length cuts (see definition in IKrawczvnski et al 
( 20061 )). and an additional cut on the arrival 



direction of the incoming gamma ray (O 2 ) re- 
ject more than 99.9% of the cosmic-ray back- 
ground while keeping 45% of the gamma rays. 
The cuts applied here are: integrated charge per 
image > 1200 digital counts («225 photoclec- 
trons), —1.2 < mean-reduced-scaled width/length 
< 0.5, and 6 2 < 0.015 deg 2 (9 2 < 0.025 deg 2 for 
the 3-telescope data set). The background in the 
source region is estimated from the same field of 
view using the "reflec ted-region" model with 8-10 
background regions (jAharonian et all 1200 lh and 
the "ring-background" model with a ring size of 
0.5° (mean radius) and a ring width of 0.175° 
(jAharonian et a]|l2005ar ). This analysis has been 
chosen to provide best sensitivity at high energies, 
in order to compare more di r ectly with the results 
reported in I Aharonian et all (|2007| ). The resulting 
energy threshold is 720 GeV. 

Observations of XMMU J063259.3+054801 (HESS 
J0632+057) were taken with the Swift satel- 
lite over four consecutive nights, from January 
26 to January 29, 2009, and were contempo- 
raneous with VERITAS observations (Set III). 
The analysis is restri cted to X-ray data from 
the XRT instrument ( Burrows et all l2005h since 
the source is not detected with BAT and the 
bright star MWC 148 causes a very high level 
of pho ton-coincidence lo sses in the UVOT instru- 
ment (jPoole et all bOOSh f 7 ! The presented Swift 



J It should be noted that MWC 148 with a B magnitude of 
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XRT observations were conducted in photon- 
counting (PC) mode with no signatures of pile-up. 
Data reduction is performed with the HEAsoft 6.5 
package. Events with grades 0-12 in the energy 
range of 0.3-10 keV are calibrated and cleaned 
using the xrtpipeline tool by applying the stan- 
dard filtering criteria and latest Swift calibration 
files. Source counts are extracted from a circular 
region with a 30-pixel radius (47.2 arcsec), and 
background events are extracted from a 40-pixel 
radius circle in a source-free region. Ancillary re- 
sponse files are generated with the xrtmkarf tool 
applying corrections for the PSF losses and CCD 
defects. The latest response matrix from the XRT 
calibration files are used in the analysis. To en- 
sure valid x 2 minimization statistics during spec- 
tral fitting, the extracted XRT energy spectra are 
rebinned to contain a minimum of 25 counts in 
each bin. 

3. Results 

Results for each of the three VERITAS data 
sets (see Table Q} , as well as for the total obser- 
vation are listed in Table 2. Figure [1] shows a sky 
map of the significance at energies above 720 GeV 
observed in the region around HESS J0632+057. 
The distribution of significances in the sky map 
is consistent with the expected distribution from 
a field with no gamma-ray source present. The 
significance at the position of HESS J0632+057 is 
2.1cr (lcr for an energy threshold of 1 TeV, see Ta- 
ble [2]) . There is therefore no significant evidence 
for gamma-ray like events from HESS J0632+057 
during the 31 hours of observations with VER- 
ITAS. The flux upper limit (E>1 TeV) for the 
complete dat a set at the 99% confidence level 
(|Helenel [l983h assuming a power-law like source 
spectrum with a spectral index of T = 2.5 is 
F(>1 TeV) < 2.6 x 10~ 13 cm^s" 1 (about 1.1% 
of the flux of the Crab Nebula; see Table [2]). This 
flux limit is ~ 2.4 ti mes lower than t he flu x re- 
ported by H.E.S.S in lAharonian et all l|2007t ): see 
Figure [2] for a light curve. The probability for a 
non- variable flux of high-energy gamma rays from 
HESS J0632+057 is derived from the VERITAS 
data and the average of the reported fluxes from 
H.E.S.S. using a x 2 -test. The test gives a \ 2 °f 



15.8 with 1 degree of freedom, corresponding to a 
probability of 0.007% (about 4<r). 

The non-detection of HESS J0632+057 by 
VERITAS initiated additional and very thorough 
quality checks of the data. Optical pointing mon- 
itors which are installed on each telescope show 
that the systematic error on the pointing was less 
than 90" during the data-taking period for Set 
II and 110. The sensitivity of VERITAS in the 
elevation range of 50-60° was confirmed by obser- 
vations of the Crab Nebula and several other weak 
gamm a-ray sources in early 2009 (e.g. lOng et al 
(|2009h ). 

The Swift XRT observations provide further ev- 
idence for X-ray flux variability in the object. The 
total photon spectrum from the four Swift XRT 
observations in January 2009 is well-described by 
an absorbed power law. Due to the low count- 
ing statistics, a fixed column density of Nh = 
3.1 x 10 21 cm ~ 2 is applied in the spectral analy- 
sis, as measured for the X MM-Newton EPIC spec- 
trum ( Hinton et a 3 l2009h . The best-fit absorbed 
power law model yields a % 2 /dof = 3.5/7, a pho- 
ton index T = 1.9 ± 0.3, and a 1 keV normal- 
ization of (2.4 ± 0.6) x 10~ 4 cm- 2 s" 1 keV" 1 . 
The deabsorbed 1-10 keV flux is (9.7 ± 2.1) x 
10~ 13 erg cm~ 2 s _1 , corresponding to a factor of 
1.8 ± 0.4 higher flux than the XMM-Newton EPIC 
data. Among the four Swift XRT observations 
no significant flux variability is measured. Fig- 
ure 3 shows the deabsorbed X-ray spectra from 
the Swift and XMM-Newton observations. Over 
the limited energy range of 0.8-4 keV, the Swift 
photon spectra is softer (r = 1.9 ± 0.3) than 
the XMM-Newton EPIC spectrum (r = 1.26 ± 
0.04). The relatively hard (T < 2) X-ray spec- 
trum is similar to the spectral slop es measured 
from the TeV b inaries LSI +61 303 (ILeahv et al 
19971 ). LS5039 (iHoffmann et all I2009TI and PSR 



B1259-63/SS 2883 (|Chernvakova et alllioOfih . No 
definite conclusion can be drawn from the X- 
ray spectral variability of HESS J0632+057. It 
appears, for exam ple from observations of PSR 
B1259-63/SS 2883 (|Chernvakova et ajl2006h that 
the photon index of the X-ray spectra in TeV 
binaries is not a strictly monotonic function of 
flux. Detailed results from the full Swift observing 



less than 9 is too faint to introduce systematic effects into 
the analysis of the VERITAS data. 



2 There were no pointing monitors installed in December 
2006 (Set I). 
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campaign of XMMU J063259.3+054801 between 
January 200 9 and April 2009 are presented in 
Falcone et al (|2009h . 

The non-detection of HESS J0632+057 by 
VERITAS provides evidence for variability in the 
flux of gamma-rays with energies above 1 TeV. 
HESS J0632+057 is unlikely to have a blazar 
counterpart. Blazars are generally hard X-ray 
point sources and bright radio sources. Neither 
the XMM observation nor the Swift images show 
the presence of any bright X-ray source other than 
the Be star within the HESS error circle. A search 
in the NVSS catalog also finds no radio source in 
the area. The VHE emission and variability can 
be easily explained if MWC 148 would be part 
of a binary system and high-energy photons are 
produced in a similar way as in LS I +61 303 
or LS 5039. The available data do not allow 
any conclusion on a possible periodicity of the 
gamma-ray signal. A detection of the compact 
companion or of the orbital motion of the Be star 
is required to confirm or refute the binary nature 
of this system. Particle acceleration and VHE 
emission from massive stars with strong magnetic 
fields has also been suggested. A confirmation 
that MWC 148//HESS J0632+057 is surrounded 
by sufficiently strong magnetic fields, along with 
further theoretical work to explain the variabil- 
ity in the gamma-ray emission, would be needed 
to establish this potentially new class of galac- 
tic gamma-ray sources. Future multiwavelength 
observations combined with results from ground- 
based and space-based gamma-ray observatories 
will provide a deeper understanding of the true 
nature of HESS J0632+057. 
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Fig. 1. — VERITAS significance map of the 
region around HESS J0632+057 for the whole 
data set and an energy threshold of 720 GeV. 
The background is estimated using the ring back- 
ground method. The location of HESS J0632+057 
is indicated by a black triangle. Also shown 
are the 95% and 99% confidence regions of the 
EGRET sources 3EGJ0634+0521 (white lines) 
and 3EGJ0631+0642 (GeV J0633+0645) (black 
lines). The Fermi source 0FGLJ0633.5+0634 is 
indicated with a 'x' sign; the blue circle d enotes 
the 95% confidence region (jAbdo et a]||2009l ). The 
circle at the bottom right indicates the angular 
resolution of the VERITAS observations. 
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Fig. 2.— The light curve above 1 TeV from HESS 
J0632+057 is shown assuming a spectral shape 
of dN/dE oc E~ r with T = 2.5. The down- 
wards pointing arrows show the 99% confidence 
limits derived here from the VER ITAS data. The 
H.E.S .S. fluxes are taken from Aharonian et all 
(|2007l ). The X-ray fluxes measured by XMM- 
Newton and Swift are indicated by open symbols. 
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Fig. 3.— X-ray spectrum of HESS J0632+057 
from XMM-Newto n EPIC data in September 2007 
(Hinto n et a 1 2009) and Swift XRT observations in 
January 2009. 
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Table 1 

Details of the VERITAS and H.E.S.S 



(|Aharonian et alI|2007D observations of 
HESS J0632+057. 



Name 


Date range 


N T ei a 


Elevation 
range 


Angular distance 
between source and 
pointing direction 


Observation 
time 
[min] 


Set I 


Dec 16 2006 - Jan 25 2007 


3 


55-65° 


0.15-0.8° 


580 


Set II 


Dec 30 2008 - Jan 03 2009 


4 


59-65° 


0.5° 


560 


Set III 


Jan 26 2009 - Jan 30 2009 


3/4 


59-65° 


0.5° 


722 


H.E.S.S. PI 
H.E.S.S. P2 


Dec 2004 

Nov 2005 - Dec 2005 








282 
372 



a Number of available telescopes 



Table 2 

Analysis results for HESS J0632+057 for E > 1 TeV. 



Upper limits <I> 7 ,ul(# > 1 TeV) are 



GIVEN AT 99% CONFIDENCE LEVEL (AFTER [HeLENfJ (Il983l )). THE INTE GRAL FLUXES AND 1 a ER RORS 
ABOVE 1 TEV REPORTED BY H.E.S.S. ARE LISTED FOR COMPARISON (lAHARONIAN ET aJ|2007| ). 





on 


off 


a 


excess 


significance 


Flux or upper flux limit 




events 


events 




events 


w 


[10- 13 cm^s" 1 ] 


Set I 


84 


594 


0.16 


-8.8 


-0.9 


< 4.2 


Set II 


131 


713 


0.16 


15.7 


1.3 


< 4.2 


Set III 


120 


669 


0.16 


11.4 


1.0 


< 3.6 


Total (I-III) 


335 


1976 


0.16 


19.4 


1.0 


< 2.6 


H.E.S.S. PI 












6.3 ± 1.8 


H.E.S.S. P2 












6.4 ± 1.5 
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